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Abstract. We demonstrate a method to determine the spectral purity of single 

photons. The technique is based on the Hong-Ou-Mandel (HOM) interference between 

a single photon state and a suitably prepared coherent field. We show that the temporal 

Q" 1 width of the HOM dip is not only related to reciprocal of the spectral width but also 

to the underlying quantum coherence. Therefore, by measuring the width of both the 

£NJ ■ HOM dip and the spectrum one can directly quantify the degree of spectral purity. 

The distinct advantage of our proposal is that it obviates the need for perfect mode 
matching, since it does not rely on the visibility of the interference. Our method is 

\/-^ , particularly useful for characterizing the purity of heralded single photon states. 
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1. Introduction 

The Hong-Ou-Mandel interference is one of the most famous phenomenon studied in 
quantum optical experiments [1]. The effect can be observed by coherently combining 
two indistinguishable photons at a beam splitter. As a result of interference in the 
Hilbert space the photons bunch together and depart at the same beam splitter output 
port. Therefore, the effect relies not only on the linear superposition of electric fields, 
but it also highlights the specific quantum characteristic of the input states. 

Originally, the HOM effect was observed by interfering two strongly correlated 
fields produced by parametric down-conversion [1-3]. High visibility was achieved by 
tight spectral and spatial filtering due to the incoherent nature of the light fields. 
Nowadays the two-photon coalescence effect constitutes the core mechanism of many 
quantum information protocols, such as linear-optics quantum computation [4] and 
teleportation [5,6]. This fact has boosted the efforts towards the experimental generation 
of pure single photon states [7], culminating in the observation of HOM interference 
between single photons emitted by a semiconductor quantum dot [8] or dissimilar 
sources [9]. 

Given the current demand on preparing high quality single photon states as resource 
for quantum information tasks, the experimental determination of their degree of purity 
and underlying spectral properties becomes paramount. At present, the standard 
adopted benchmark is the visibility of the HOM interference. However, since it indicates 
the degree of indistinguishability between the quantum states [10, 11], it renders a 
pessimistic estimate for the photon purity as it is highly sensitive to mode mismatching. 
Although a complete mapping of the spectral properties could be realized [12,13], we 
are here concerned with the characterization of a single quantifier, i.e. the degree of 
purity. 

As already noticed in the seminal Hong-Ou-Mandel paper [1], the fourth-order 
interference experiment provides a way to determine properties of the interfering photons 
with extremely high temporal resolution, on the order of femtoseconds. In the earlier 
experiments, the temporal width of the dip was mainly determined by the inverse 
bandwidth of the employed spectral filters, which was indeed equal to the coherence 
time. This fact was just a consequence of the choice of input state-an incoherent 
parametric down-converted one. In a more general situation the width of the HOM 
dip can greatly differ from the coherence time. Still, the properties of the temporal 
width of the HOM interference have been mostly ignored. We show here that the latter 
is related to both the reciprocal of the spectral width and the underlying quantum 
coherence, thus providing crucial knowledge about the state's spectral purity. 

In this paper we present a method to characterize the degree of spectral purity in 
a robust way and with few measurements. In order to do that we analyze the temporal 
profile of the HOM interference between a single photon state and an auxiliary weak 
coherent field, used as a reference. The advantage of our proposal is that it enables the 
quantification of the purity without the need for perfect mode matching, thus greatly 



Accessing the purity of a single photon by the width of the HOM interference 3 

reducing the technical demands on the experimental conditions. In the simplest case, 
our method requires only the knowledge of the spectral widths of the reference and single 
photon state, together with the temporal width of the HOM dip. The first two quantities 
are in general easily measured by standard high-sensitive spectrographs [14,15]. 

Our technique gives exact results when two assumptions commonly found in the 
literature are made. Moreover, it provides good results even in the most general case. 
Our proposal is well adapted for photon pair sources, which we discuss theoretically, 
showing the expected performance under different conditions. The best situation occurs 
when employing a practical and popular strategy to generate pure states, the one relying 
on tight spectral filtering. For other situations an error limited to about 10% may occur 
in a worse case scenario. To provide concrete results that illustrate the applicability of 
our technique, we also present an experimental investigation. Our findings corroborate 
the usefulness and benefits of the method. 

This paper is organized as follows. In section 2 we present general results relating 
the width of the HOM dip with the spectral purity of the single photon state. We define 
the photon spectral density function and show that in a Gaussian scenario it can be 
parametrized by two widths. While the major one is the spectral width, the minor one 
is related to the reciprocal of the HOM time duration. The purity acquires a simple 
geometrical interpretation, and is shown to be the ratio between those two widths. 
The first order coherence time, on the other hand, is the inverse of their difference. 
Therefore, as well known, spectrally pure states have a coherence time that tends to 
infinity. Impurity, by contrast, leads to a broader spectrum, thus shorter coherence 
time, resulting in a wavepacket that is not Fourier-transform-limited. After establishing 
those concepts we introduce our method for characterizing the purity. In section 3 we 
discuss the validity of our proposal when applied to states arising from photon pair 
sources. Section 4 covers an experimental investigation, where we employ our technique 
to characterize the spectral purity of signal and idler photons generated in a wave- 
guided parametric down-conversion source. Further theoretical details about photon 
pair sources is provided in the Appendix. 

2. Purity, Temporal Coherence and HOM Interference 

Before explaining how to retrieve the spectral purity from the HOM interference, we 
first define a pure single photon state prepared with a spectral shape u(u) as [16] 

|tf) = [ duu*(uj)P(u)\0). (1) 



From its corresponding density operator one observes that the spectral density function 
is simply g pme (u>, u/) = u* (u>)u(u>') , which is symmetric with respect to u> and u/. 
Considering a field with a Gaussian spectral shape the function g pure presents a circular 
symmetry in frequency space, as illustrated in figure 1 (left). To include the description 
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Figure 1. Contour lines of the density function \g(u, o/)|. Left: Pure state with 
Gaussian spectral profile. Right: Pictorial representation of an impure spectral density 
function, with major and minor widths denoted by o g \ and a g 2, respectively. This 
particular form of g is separable along the ±45° axis. 



of impure states, a more general representation for the single photon state is used 

p s = du du' g{u,u') a j (u)\0)(0\a{u') . (2) 

Since density operators are Hermitian, any spectral density function g is symmetric 
about the 45° axis, i.e. about (oj + to'). As usual, the diagonal elements of the density 
matrix correspond to the populations. In the present case those compose the spectrum 
of the field, which is directly measured with a spectrograph and whose width we denote 
by cr g i, as depicted in figure 1 (right). As a final remark, we recall that the degree of 
purity is calculated as 



V s = Tr(p^) = du duj' \g(u,u') 



(3) 



Therefore, the shape of \g(u,u')\ offers direct intuitive insight into the purity of the 
quantum state. In particular, if \g(u,u')\ is a 2-dimensional Gaussian function the 
purity is just the ratio between minor an major axis, as derived later. 

To provide further intuition on those physical concepts, we can also calculate the 
first order correlation function [17] 



G {1) (tu,co') = Ti[p s dl(co)d(uj')} = g(u,u') 



(4) 



We realize that the spectral density function is equal to the non-normalized first order 
frequency correlation function. The condition of purity as the separability of g in the 
{u,uj'} frequency space corresponds to the condition of perfect first order coherence 
as defined by Glauber [17]. Thus, the techniques for tailoring the spectral properties 
of single photons also apply for engineering wavepackets with a specific temporal 
coherence [18]. 

The method we propose for characterizing the spectral purity V s relies on the Hong- 
Ou-Mandel interference between a single photon state p s (called signal) and a highly 
attenuated coherent state |/3) with a spectral shape u{u) (called reference) [19-21]. 
The experiment is schematically shown in figure 2. Both fields interfere in a 50/50 
beam-splitter, after which avalanche photo diodes (APDs) record the number of click 
events. The measurement of coincident events versus the time delay r between signal and 
reference fields reveals the HOM interference dip. The reference field can be provided 



Accessing the purity of a single photon by the width of the HOM interference 5 

1/3} 






APDi 



EHT 



APD 2 



Coincidence 
discriminator 



Figure 2. Scheme for realizing the Hong-Ou-Mandel interference between two fields 
with a relative optical delay r. APD: Avalanche photodetector 



by any independent laser, since there is no need for phase stabilization in this kind of 
interference. 

Given that the coherent state \0) contains higher photon number contributions, to 
calculate the probability of coincidences we consider the weak power regime in which 
the interference between the reference and signal comes mostly from single photons 
in each mode. So far we have ignored the photon number purity of the signal state. 
Generally, this one is well characterized by the second order coherence [22,23], which 
can be measured with a Hanbury-Brown-Twiss interferometer. Particularly for photon 
pair sources, it is known that the prepared state presents a small probability of two 
photons. We include the degrading effect of multiphoton components in a first order 
approximation, i.e. they are assumed not to interfere resulting only in background 
events. 

The joint probability of registering photons at detectors APDx and APD 2 at times 
t\ and t 2 is given by 

Pia(ti,ta) = (Et\t 1 )Et\t 2 )Ei + \t 2 )E[ + \t 1 )) , (5) 

where the average is calculated over the state of the single photon source and the 
reference field. The fields E\ and E 2 are related to the reference Ei and signal 
Eg by the beam splitter transformation 

E^ih) = -= / dco 1 [a s (co 1 )e~'^ r + a^^e^^, (6) 

V 2 Jo 

E { 2 + \t 2 ) = — / du 2 [d s (u 2 )e-^-ap(u 2 )}e-^. (7) 

V 2 Jo 

Since the detection response time is much slower than the correlation time between 
the fields, we integrate Pi 2 (ti,t 2 ) over all possible times. This leads to a probability of 
coincidences given by 

P c (r) = B + S[l - l(r)} , (8) 

with 

X{t)= [ [dco dJ u*{u) g(co, to') u{u')e iT{w - ul) . (9) 



The parameter B is a constant proportional to the probability of vacuum and two 
photons in both signal and reference fields, whereas S is proportional to the single 
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photon contributions [24] . The function X(t) accounts for the interference between the 
single photon components of both fields. Its value at zero delay T = 1(0) gives the 
overlap between the spectral density functions of signal and reference. 

The visibility of the HOM interference is defined as the fractional reduction of the 
probability of coincidence from its uncorrelated value, i.e. 

P c (r ->■ oo) B + S { ] 

We note in particular that, by removing background events from P c (r) one obtains a 
"corrected" visibility that is equal to the spectral overlap factor T. Experimentally, 
mode mismatch in other degrees of freedom decrease even further the value of V from 
the best one T. We emphasize that the latter is only constrained by the spectral overlap 
between the density functions, which includes the property of mixedness in the signal 
state. 

One observes from (9) that perfect interference is only achieved when the spectral 
density function g(uu,u>') is separable, i.e. the signal photon is pure, and matches the 
spectrum of the reference u(u). In this special case, X(r) is nothing but the convolution 
of the temporal distributions of reference and signal. In a general situation, however, 
the width of the HOM dip will depend not only on the diagonal elements of g(u, a/) but 
on the complete spectral density function. It is precisely the shape dependence of the 
latter with the purity that we explore. 

2.1. Method for Characterizing the Spectral Purity 

Our aim is to find a simple relation between (3) and (9). Given the symmetry of this 
problem, we analyze it in a 45° rotated frequency basis and denote by {x, y} the new 
reference frame. The latter follows the direction of the major and minor axis of g(u, u/), 
as illustrated in figure 1 (right). Our first assumption consists in considering that the 
spectral density function of the analyzed state is separable in the {x, y} basis, thus 
described by g(u,u f ) = gi(x)g 2 (y). In this case the purity and the HOM interference 
dip are respectively given by 



V s = j j dxdy\ gi (x)\ 2 \g 2 (y)\\ (11) 

X(r) = [ dxu*(x) 9l (x) [dyu(y)g 2 (y)e tT y. (12) 



We observe that the temporal shape of X(r) is only affected by the elements of the 
spectral density function lying at -45°. Performing the Fourier transform of the HOM 
dip allows one to recover the function g 2 (y). In particular, if both X(r) and the spectrum 
of the reference are Gaussian functions, then g 2 (y) must also be. We recall that the 
spectrum of the studied signal is determined by the elements of the g function lying at 
+45°, in other words by gi(x). Therefore, the shape of the latter can be verified with 
an spectrograph. 
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Our second assumption consists in approximating the spectral density function by 
a Gaussian distribution and considering only linear phase terms. This approximation 
is valid in many experimental situations, as we discuss later on. Other possibility, 
more realistic for quantum dot sources would be to assume a Lorentzian profile. Such 
modification would lead to similar conclusions. Under this condition, 

1 f x 2 y 2 



g(x, y) = exp -— -^ - — -^ + iny ) , (13) 

'2nal V 4(T 9 i 4cT 9 2 



with k a constant that depends on the physical system under consideration. For a 
photon pair source the phase term depends only on the difference of frequencies "y". 
The normalization factor takes into account the condition of positivity that any physical 
density operator has to satisfy. With this definition of g, the standard deviation of the 
spectrum (field amplitude) is precisely a g i, and the spectral purity is given by 

V a = o- g2 /o- g i ■ (14) 

Lack of purity acquires a simple geometrical interpretation: the g function becomes 
compressed, displaying an elliptical shape in the frequency space. To characterize the 
purity one needs to access the values of a g \ and a g 2- While we have already discussed 
that the first can be measured with a spectrograph, for the second one we analyze the 
HOM interference. 

To solve (12), we consider a reference field with a Gaussian spectrum characterized 
by a standard deviation o~$. The HOM dip is determined by X(r) = T exp [— r 2 /(2<5 2 )], 
with the following overlap factor 



T=Tihr- 1> ( 15 ) 



2 2 ' 



and temporal width 



^ + ^J- (16) 

First, we emphasize that the overlap factor characterized via the width of the HOM dip 
quantifies exclusively the amount of spectral overlap between u*(lo)u(lo') and g(u>,u>'). 
This would not be the case when employing the visibility via equation (10), given its 
sensitivity to mode mismatch in all degrees of freedom. Second, the width of the dip 
enables the characterization of a g2 - Since a g2 = o~ g \P s the width of the HOM dip is not 
only related to reciprocal of the spectral width, as generally believed, but also to the 
purity of the single photon quantum state. Hence, the knowledge of very few-directly 
measurable-parameters completely determines the purity of the interfering photon, 

V s = -. l ^=- (17) 

y/2o* gl [S 2 - l/(2o%)) 
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For completeness, we also calculate the coherence time t c of the single photon state. 
Using the generalized Wiener-Khintchine theorem [25], the temporal correlation function 



is determined by -f(t,t') = r(t,« / )/v r (*>*) r (* / >*')> with 

T(t, t')= f f du dio'g(co, u'y^-"'^ . (18) 

The coherence time is proportional to the standard deviation of ^y(t,t'), which results in 
t c oc l/(o"pi — o~ 2 g2 ). Thus a pure state has a coherence time that tends to infinity, while 
incoherent chaotic sources have a coherence time that is equal to the inverse of their 
spectral bandwidth. 

3. Application to Photon Pair Sources 

By employing nonlinear processes such as three- and four-wave mixing one can generate 
photon pairs from an intense pump beam. The fundamental physical constraints 
governing these processes are energy conservation and phase- matching of pump and 
down-converted (signal and idler) waves. The first results in a strict photon-number 
correlation between the twin photons. Therefore, one of them can be used to herald 
the existence of the other. Moreover, those conservation laws lead to a two-photon 
state characterized by a complex, highly correlated joint-spectral function [26]. In our 
context, the spectral correlations together with the fact that the triggering photon is 
measured with a conventional frequency-blind APD, lead to an effect akin to dephasing 
and the heralded photon ends up in a mixed state. One possibility to produce highly 
pure single photon states relies on careful engineering of the phase-matching condition 
of the non-linear material [27]. Progress in this direction has been described in several 
recent experiments [28-30]. Nevertheless, a more common strategy to generate pure 
states still relies on spectral filtering. 

We analyzed the performance of our method considering a single photon state 
produced by parametric down-conversion. To follow our analysis it is only necessary to 
know that the joint-spectral function <f>(u s , Vj) of the photon pair state is the product of 
the phase-matching (PM) condition $(z/ s ,z/j) and the pump envelope a(u s + Vi). The 
first has a spectral shape determined by a sine function and a slope 9, the latter is 
assumed to be a Gaussian. The parameter v s (z/j) denotes the frequency detuning of 
signal (idler) with respect to the central frequency of the pump. Figure 3 allows to 
quickly overview the relevant parameters. We recall that the spectrum of one field 
is just the marginal distribution obtained by integrating <fi over the frequencies of the 
other. Finally, considering that the idler photons are used as a trigger, the spectral 
density function of signal is given by [24] 



U s , l/) OC / dVi (j)(l/ 3 , Vi)<f>*(l/ a , Vi)\/ faiysW fail's) fifa) , (19) 



9[ 

with f^u^ the intensity transmission function of the applied spectral filter. We note 
that the joint-spectral function <fi determines the frequency correlations between signal 
and idler, being a property of two fields. The density function g, although connected to 
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Figure 3. Illustration of the joint-spectral distribution (ellipse), determined from the 
properties of the pump beam (a) and phase-matching ($). The standard deviation of 
the pump spectral function and PM are respectively denoted by a p and tr$. In what 
follows we will use the notation w to refer to the full width at half maximum (FWHM) 
of pump (p), phase- matching (<fi), spectral filter (/) and reference field (/3), which is 
indicated by the appropriate subscript. 



4>, is a property of a single beam and has a meaning directly connected to the degree of 
spectral coherence. Thus care must be taken to avoid confusing these two functions. 

To test our technique, we simulated several conditions, resulting from different 
choices of the joint-spectral function and eventual use of Gaussian spectral filters. Those 
define the spectrum and purity of the signal field. Taking a reference beam with similar 
spectral width we evaluate (9) for several temporal delays, resulting in a specific HOM 
dip. Fitting Gaussian functions to the fields' spectra and HOM dip, we obtain the 
necessary parameters to evaluate the purity via (17), which will be denominated as 
"Pwidth- Finally, we explicitly compare the value obtained from our method with the 
one found from a direct calculation of the purity as V=Ti{p 2 ] (3). In what follows we 
considered the HOM interference of signal photons, thus g will never refer to the density 
function of idler. 



3.1. Performance of the method with different wavepackets 

As a general rule, the shape of the HOM dip is dominated by the function with the 
smallest spectral width. Thus a very narrow phase-matching (PM) function leads to 
a squared profile for the temporal dip (Fourier transform of the sine). If on the other 
hand the width of the pump field is decreased such that it becomes on the order or 
narrower than the PM (in frequency units), then the shape of the HOM dip tends to 
a Gaussian. Below we discuss the main factors to be considered when applying our 
method, illustrating our proposal with concrete results for three different situations of 
great interest. 

(I) Highly frequency correlated joint-spectrum: In this case, to increase the degree 
of spectral purity of the heralded single photon it is necessary to apply tight filtering. 
High purity is achieved when the filter bandwidth is on the same order as the intrinsic 
width of the joint-spectrum (minor axis). We show in figure 4 (left) a filtered joint- 
spectral function that closely resembles the one found in our experiment, which will be 
discussed later. The corresponding spectral density function of the signal \g(u, u')\ is 



Accessing the purity of a single photon by the width of the HOM interference 

3a 1 



10 



> 






* 



-2 2 
Signal Wavelength (nm) 





-2 2 

Delay (ps) 



Figure 4. (Left) Filtered highly correlated joint-spectral function. (Middle) Contour 
lines of \g(v, v')\ as a function of wavelengths. (Right) HOM dip obtained by interfering 
the signal with the reference field: (circles) calculated using (9), (solid line) Gaussian 
fit. The following parameters were used: w p = 2 nm, w$ — 0.8 nm, 9 = 55°, Wf = 1 nm 
and wp = 1 nm. 



presented in figure 4 (middle). Its shape is very close to a function that is separable 
along the major and minor axis. Employing the width of the HOM dip, shown in figure 4 
(right), one finds the degree of purity estimated with our method: Pwidth = 59%. This 
value has to be compared with the spectral purity obtained by the direct calculation: 
V = 60%. We find good agreement. In general, to increase the purity of the signal 
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Figure 5. (Left) Width of the HOM dip (defined by the FWHM of a Gaussian fit) 
for different widths of the phase-matching, considering unfiltered and 1 nm filtered 
triggering photons. (Right) Comparison between "true" purity (solid line) and the one 
given by our method (dashed line) for unfiltered (dark colours) and 1 nm filtered (light 
colours) triggering photons. The following parameters were used: w p = 2 nm, 6 = 55°, 
Wf — f nm and wp = 1 nm. 



state even further, a second filter can be included in the idler field. This non-local 
manipulation of signal's frequency leads to a narrower heralded spectrum, but at the 
same time to a narrower HOM dip. The latter is a consequence of describing the phase- 
matching by a sine function and would not be observed in case it could be approximated 
by a Gaussian. This effect can be seen in figure 5 (left), where we present the width of 
the dip for different PM widths. Thus, by adding a 1 nm filter to limit idler frequencies 
one finds V = 89% and "P w idth = 88%. Again our method proves to be a faithful tool for 
the characterization of the purity. Its efficacy as a function of the PM width is presented 
in figure 5 (right). For small values of the PM width the resulting squared HOM dip 
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profile indicates that the #2 function is not a Gaussian. This leads to retrieving a value 
for the purity that is about 10% underestimated. As a general behaviour, the estimate 
"Pwidth tends from below towards the correct value of the purity, being optimal when 

W$ « Wf. 

(II) Phase-matching lying at 90°: High purity heralded single photons are generated 
by appropriate choice of the pump width and there is no need for spectral filtering. A 
typical joint-spectrum and the corresponding density function of the heralded photon are 
illustrated in figure 6 (left) and (middle), respectively. In this case, by tracing over idler 
frequencies to obtain p s one naturally introduces a smearing out of the sine structure, 
resulting in a g(is, v') function that is a 2-dimensional Gaussian. Therefore, our method 
can be used without restrictions and we find that V = "Pwidth = 96%. 
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Figure 6. (Left) Uncorrelated joint-spectral function. (Middle) Contour lines of 
\g(is, v')\ as a function of wavelengths. (Right) HOM dip obtained by interfering signal 
with the reference field: (circles) calculated using (9), (solid line) Gaussian fit. The 
following parameters were used: w p = 1.8 nm, w§ — 1 nm, 9 — 90° and wp = 10 nm. 

(Ill) Phase-matching lying at 45°: Again high purity is achieved by matching the 
widths of pump and the PM, resulting in a condition that is the most unfavorable for our 
technique. The reason is that the "wings" of the sine function appear in a pronounced 
and symmetric way, thus surviving to the smearing effect discussed above. We show in 
figure 7(left) and (middle) a typical joint-spectrum and density function, respectively. 
Although the function g(y, is') is not separable and the validity of the method breaks, 



© 



-5 5 

Signal Wavelength (nm) 




-4-2024 
A (nm) 




-0.5 0.5 

Delay (ps) 



Figure 7. (Left) Uncorrelated joint-spectral function. (Middle) Contour lines of 
\g(is, v')\ as a function of wavelengths. (Right) HOM dip obtained by interfering signal 
with the reference field: (circles) calculated using (9), (solid line) Gaussian fit. The 
following parameters were used: w p = 1 nm, w$ = 3 nm, — —45° and wp — 4 nm. 



a blind application of (17) results in a characterization of the purity that is only about 
12% wrong, nevertheless in this case it is overestimated. This state has a purity of 



Accessing the purity of a single photon by the width of the HOM interference 



12 



V = 88% and, given its reduced value, the use of traditional spectral niters is still 
inevitable. The local application of a Gaussian spectral filter to the single photon state 
effectively smooths the effect of the sine bringing the spectral density function towards 
the ideal conditions of our method. Applying a 5 nm filter to the signal field, enough 
for cutting most of the "wings" arising from the sine function, would result in V = 96% 
and a 4% overestimated P w ;dth- 

4. Experiment 

To show the experimental feasibility of the method, we measured the HOM dip between 
one of the twins and a coherent reference. An illustration of the experimental setup is 
presented in figure 8. A 1.45 mm long, type-II, periodically poled KTiOP0 4 waveguide 
is pumped by a frequency doubled Ti: Sapphire laser operating in the ultrafast regime 
(796 nm, FWHM equals to 10 nm, autocorrelation length of 170 fs and repetition rate 
of 4 MHz). A small fraction of the fundamental light is used as the reference field. The 
spectral width of the pump field is w p = 2.0(1) nm. The path delay c r between reference 
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Figure 8. Sketch of the experimental setup. SHG: second harmonic generation, SF: 
spectral filter (wf = 1.0(1) nm), HWP: half-wave plate, PBS: polarizing beam-splitter, 
APD: Avalanche photodetector. 



and one of the twins is adjusted by using a high-precision translation stage, which is 
moved in steps corresponding to a delay of 0.22(1) ps. The fields are combined by using 
a fibre integrated beam-splitter, thus ensuring a good spatial overlap. We measure 
coincident clicks using two avalanche photodetectors whose signals are acquired by a 
computer using a time-to-digital converter (TDC). 

In our case the joint-spectrum is highly correlated and we prepare heralded single 
photons by applying 1.0(1) nm spectral filters to both, signal and idler photons. 
Discarding any information about idler (unconditional measurement), it is legitimate 
to assume that the spectral width of the signal is equal to the filter width o g \ = a/, 
since the unfiltered marginal spectrum is much broader. A similar condition holds for 
the reference beam, cjp = af. Thus, to avoid measuring a g i, we begin the experiment by 
investigating the interference between the unconditional signal state and the reference, 
i.e. between a thermal and coherent state. 

The result of our measurement is presented in figure 9, where we show typical HOM 
dips observed by interfering either signal or idler with the reference beam. By analyzing 
several similar sets of data we extract the following FWHM values: w s = 1.75(5) ps and 
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itfj = 1.54(5) ps. The purity of the unconditional mode is retrieved applying (17). One 
finds that P„ idth = 0.64(10) for the signal and, in a similar way, "P„/ dth = 0.77(10) for 
the idler, despite the low visibilities (depicted in figure 9). 
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V = 21% 
FWHM= 1.72(7) ps 
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-3 -2 



-1 1 
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Figure 9. Measured HOM dip by interfering the reference field with: (Left) signal, 
(Right) idler. The width of the dip and the corresponding visibility are shown in the 
plots. 



P0 (%) Pi (%) V2 (%) 



Signal 99.7015 0.2978 0.0006 0.0030 
Idler 99.8309 0.1689 0.0002 0.0027 



Table 1. Mean photon number of the reference field and photon number statistics 
of the field that interferes with it. 



The best way to compare our method with the standard one based on the visibility 
is to analyze the value of the overlap factor, since for the interference between two 
intrinsically different sources the visibility is not directly connected to the purity [10]. 
As discussed, this quantity can be retrieved either from the HOM width via (15) or 
from the visibility via (10). Here we refer to the former as T and to the latter as 
T vis . For the first case we directly obtain T^ s > = 0.76 and T® = 0.86. For second 
case, we measured the photon number distribution of both fields using a standard 
time- multiplexing detector [24]. In our experiment we can neglect the effect of dark 
counts by applying a tight temporal gating to the acquisition (~3 ns). The observed 
probabilities of vacuum p , one-photon p\ and two-photon p 2 components of the studied 
state, together with the mean photon number |/3| 2 of the coherent field, are shown in 
Table 1. With these values one can calculate S = p\ |/3| 2 and B = [po |/3| 4 /2 + p 2 ], which 
yield an overlap factor of T v f s = 0.46 and T vl * s = 0.60, smaller than the ones retrieved 
from the width of the dip. This exemplifies well the fact that the visibility is affected by 
mode mismatch and other experimental imperfections, contrarily to the temporal width 
used in our method. 

Concerning the heralded photon state, the non-local spectral filter in the idler arm 
reduces the spectral width of the signal state o~ g i, thus differing from o~f (bandwidth of 
the filter in the signal arm). Therefore, in this scenario a g \ has to be measured in order 
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to allow the characterization of the heralded purity via the HOM width. Unfortunately, 
in our case this width is narrower than 1 nm and can not be directly accessed in our 
system. Still, the degree of purity can be estimated making use of the knowledge that 
our source produces photon pairs, described by a joint-spectral function, with measured 
unconditional purities. In previous measurements we characterized the slope of the 
phase-matching function to be 9 = 54.7° [24], thus by comparing the values of V^ dth and 
^width with a theoretical curve (obtained by numerical integration), we can characterize 
the PM width. From this analysis we estimate that, employing a 1 nm filter in the idler 
arm, the heralded purity is Vh = 0.90(6) and the spectral overlap with respect to the 
same 1 nm broad reference field is T = 0.88. 

5. Conclusion 

We demonstrated a new method to directly determine the purity of single photons in 
a way that is independent from mode matching. We employ the standard Hong-Ou- 
Mandel interference between two fields but, contrarily to the usual approach based 
on the amplitude of the interference, we benefit from the robustness of its temporal 
width against mode mismatch. Its applicability for characterizing the state produced 
by photon pair sources is explicitly examined and discussed for different experimental 
scenarios. In particular, we show its advantage in an experiment dealing with tightly 
filtered parametric down converted photons. 

Our technique can also be extended to situations in which signal and reference 
with the same wavelength are not easily available. In this case, the solution might be 
to employ two independent but equivalent photon pair sources and interfere the signal 
mode from each of them. The information about the purity should still be available in 
the temporal width of the HOM dip. 

We hope that our work will trigger further investigations about more practical and 
robust techniques for quantifying the quality of general single photon states. We believe 
that our ideas will boost the characterization and optimization of high-purity heralded 
single photon states, thus having direct impact on the practical realization of quantum 
information protocols [31-33]. 

Appendix 

The precise expressions used in our simulations are provided below. The phase-matching 
is defined as 

$(u„ Ui) = smc(Ak(co s , u i )L/2)e lAk(u)s ^ )L/2 , (20) 

where {u s ,Ui} represents the frequency space of signal and idler fields, Ak the phase- 
mismatch between pump and twins and L the length of the nonlinear medium. We 
define v^ (/j, = s,i,p) as the detuning relative to the central frequencies of signal, idler, 
and pump, respectively. Moreover, the first order expansion of the phase-mismatch 
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reveals that Ak(v s , z/j) = k s v s + Kji/j, where k^ is the respective group velocity mismatch 
with respect to the pump. With this consideration the PM function is rewritten as 

*(!/., v % ) « sine (<**0v.+^0A e «#(^. + ^) , (21) 

with tan# = K s /Ki, cr^ 2 = (k 2 + K 2 )'jL 2 /2, and 7 = 0.193. Since directly measurable 
observables are proportional to the intensity, we define the phase-matching width w$ 
by the FWHM of |$| 2 , i.e. w$ = 2(j $ >/Ei2. 

The pump field is considered to have a Gaussian spectral profile, defined by 

a{u 9 , Ui) = exp(-(i/, + i/i) 2 /(2<Tp)) , (22) 



with FWHM (intensity) given by w p = 2cr p vm2. Finally, using (21) and (19) we note 
that the linear phase term of g(zA;, 1/) is [L/2 « s (za, — i/ s )]. 

In part of our analysis we employed spectral filters. Those are defined by Gaussian 
functions f(u), such that f{oo) 2 has standard deviation and FWHM respectively denoted 
by o~f and Wf. Moreover, for the characterization of the purity we assumed a reference 
field with Gaussian spectral profile u(u). The standard deviation of u(u) 2 is denoted 
by o~p and the respective FWHM by wp. 
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